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ABSTRACT 

Optical and near-infrared (NIR) spectra are presented for ULAS J13191 1.29+095051.4 (hereafter ULAS J1319+0950), a new redshift 
z = 6. 127 ± 0.004 quasar discovered in the Third Data Release (DR3) of the UKIRT Infrared Deep Sky Survey (UKIDSS). The source 
has Fveea = 19.10 ± 0.03, corresponding to Mi45o,ab = -27.12, which is comparable to the absolute magnitudes of the z - 6 quasars 
discovered in the Sloan Digital Sky Survey (SDSS). ULAS J1319+0950 was, in fact, registered by SDSS as a faint source with 
Zab = 20. 13 ± 0. 12, just below the signal-to-noise ratio limit of the SDSS high-redshift quasar survey. The faint z-band magnitude is 
a consequence of the weak Ly a/N v emission line, which has a rest-frame equivalent width of ~ 20 A and provides only a small boost 
to the z-band flux. Nevertheless, there is no evidence of a significant new population of high-redshift quasars with weak emission lines 
from this UKIDSS-based search. The Lya optical depth to ULAS J1319+0950 is consistent with that measured towards similarly 
distant SDSS quasars, implying that results from optical- and NIR-selected quasars may be combined in studies of cosmological 
reionization. 

Also presented is a new NIR-spectrum of the previously discovered UKIDSS quasar ULAS J020332. 38-1-001229. 2, which reveals the 
object to be a broad absorption line quasar. The new spectrum shows that the emission line initially identified as Ly a is actually N v, 
leading to a revised redshift of z = 5.72, rather than z = 5.86 as previously estimated. 

Key words, quasars: individual: ULAS J020332.38+001229.2; ULAS J13191 1.29+095051.4 - infrared: general - cosmology: ob- 
servations 

1. Introduction matic discovery from these studies is the marked increase in 

I I J . I I the optical depth to neutral hydrogen at redshifts of z ^ 5.7 

Since their discovery by |Schmidl| (|1963[) and [Hazard et alj (Becker et al. 2001; Fan et al. 2002, 2006a). The increase in op- 
1119631), quasars have continued to be t he most rev ealing probes jjcal dept h appears to represent th e end of cosmological reioniza- 
of the high-redshift Universe (e.g., ISchneideij [199£). Even tion (e.g.. lBarkana & Loebl2001h . a conclusion su pported by the 
though galaxies have been detected out to greater di stances (e.g., Wilki nson Micmwave Anisotropy Probe (WMAP; iBennett et al.1 
IOK-1, with a redshift of z = 6.96. 'Ive et al]l2006h and gamma 12003') measurements of the cosmic microwave background (e.g., 
rai^bursts may briefly be more luminous (e.g., [Haislip et al.| Punklev et al. 2009). While a consistent picture of reionization 
|2006|), the highest-redshift quasars (e.g., SDSS 1 148+525 1 at has emerged, the direct measurements of this process are Umited 
z = 6.42 , iFan et alj|2003t CFHQS J2329-0301 at z = 6.43, to the smafl number of z 6 quasars known, and it is clear that 
IWillott et al. 2007) are more useful because they remain bright jhe discovery of quasars with z ^ 7 is vital to further progress in 
enough to be investigated in detail. It has been possible to ^jjis field 
obtain high signal-to-noise ratio (S IN) spectra for all the 

known z - 6 quasars, robustly confirming the na ture of these The majority of the known z - 6 quasars have been identi- 
sources, revealing their intr insic properties (e.g.. IWalter et alj fiedby looking for point-sources with veiy red optical colours in 
120041; IVenemans et alj|2007l) and, through absorption, probing wide-field survey s such as the Sloan Digital Sky Survey (SDSS; 
the intervening matter out to the quasars' redshifts. The resul- lYork et a T 2000*) and the Canada France High-z Quasar Survey 
tant measurements of the z - 6 quasar population are also of (CFHQS; Willott et al. 2007), and more discoveries will be made 
interest, as they provid e critical hmits on the early structure for- in such projects. Optical searches are, however, unlikely to probe 
mation scenarios (e.g.. lKurk et al.ll2007i) . Probably the most dra- beyond the cuiTent redshift limits. Almost all z - 6 quasar emis- 
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sion at rest-frame wavelengths shorter than the Ly a transition 
at A - 0.1216/im is absorbed by intervening hydrogen, and 
such quasars are effectively dark at observed wavelengths below 
A ^ [0.85 + 0.12(z - 6)];um. Conversely, most optical charge- 
coupled device (CCD) detectors have a poor response beyond 
wavelengths of A ^ I yum (i.e., redward of the z or Z bands), so 
quasars with a redshift of z ^ 6.4 are destined to remain invisible 
to CCD-based surveys. Given the rarity of z ^ 6 quasars (e.g., a 
surface density of ~0.02deg"^ to zab = 21.0: lJiang et al.ll2007l) . 
progress can only be made using wide-field surveys at longer 
wavelengths. 

In the long term, radio surveys with, e.g., the Low Frequency 
Array (LOFARQ) and the Square Kilometre AiTay (SKAB will 
prov ide a powerful complementary approach (e.g.. lWvithe et al.l 
12009 ). but the first steps beyond the current limits will come in 
the near-infrared (NIR), with surveys following the same basic 
principles as the SDSS and CFHQS searches. The largest com- 
pleted NIR survey, th e Two Micron All Sky Survey (2MASS; 
ISkrutskie et al.ll2006h . with a magnitude limit of 7vega - 15.8, 
does not have suflicient depth to find any plausible high- 
redshift quasars. The Vis ible and Infrared S urvey Telescope 
for Astronomy (VISTA; lEmerson et al.l |2004|) should cover ~ 
2x10"* deg^ to 7vega - 20 during the next decade, but progress 
in the search for high-redshift quasars will come first from the 
partially complete UKIRT Infrared Deep Sky Survey (UKIDSS; 
[Lawrence et al. 2007). 

One new high-redshift quasar, 

ULAS J020332.38-I-001229.2, hereafter UL AS J0203-H0012 
with an estimated redshift of z = 5.86 jVenemans et al.l 
|2007|) . has already been discovered in UKIDSS; the sec- 
ond such discovery, ULAS J1319n.29-t-09505L4, hereafter 
ULAS J1319+0950, is presented here. Section |2] gives an intro- 
duction to the UKIDSS LAS project and Section |3] describes 
the selection techniques that led to ULAS J13 19+0950 being 
identified as a candidate high-redshift quasar Optical and NIR 
spectra of ULAS J1319+0950 are presented and compared 
to those of z - 6 SDSS quasars in Section |4] In addition, 
a new NIR spectrum of ULAS J0203 -1-0012 is presented in 
Section |5] together with a revised redshift estimate for this 
source. The conclusions and future prospects for high-redshift 
quasar searches with UKIDSS are discussed in Section|6] 

All photometry is given in the native system of the tele- 
scope in question, and explicitly subscripted. Thus SDSS / and 
z photometry is on the AB system, whereas UKIDSS Y and J 
photometry is Vega-based. The AB corrections for the UKIDSS 
bands are Fab - i'vega + 0.634 and 7ab - Jvega + 0.938 
dHewett et al.ll2006l) . Calculations of absolute (AB) magnitudes 
are performed assuming a fiducial flat cosmological model with 
normalised matter density Qm = 0.27, normalised vacuum den- 
sity Qa = 0.73, and Hubble constant Ho = 71 kms"' Mpc"'. 

2. The UKIRT Infrared Deep Sky Survey 

UKIDSS (iLawrence et al ] l2007h is a suite of five surveys under- 
taken with the Wide Field Camera (WFCAM; ICasali et al.l2007h 
on the 3.8 m United Kingdom Infrared Telescope (UKIRT) at 
Mauna Kea, Hawaii. The WFCAM detectors are four sparse- 
packed 2048 X 2048 Rockwell Hawaii-II arrays, each of which 
has a field of view of 0.05 deg^. Observations are generally un- 
dertaken in sets of four contiguous pointings which, together, 
cover 0.77 deg^. Full details of the survey operations can be 

' See the LOFAR web-site at http : //www . lof ar . org/. 

^ See the SKA web-site at http : //www . skatelescope . org/. 



found in iDve et al.l (l2006h . The individual i mages are analysed 
by the data reduction pipeline described by 'Ir win et all (l2009l) 
and the catalogues of detected objects are merged across bands 
into a queryabl e relational database at the WFCAM Science 
Archiv^ (WS A; iHamblv et al.ll2008h . 

The UKIDSS Large Area Survey (LAS) will cover 4000 deg^ 
within the SDSS footprint with a series of 40 s exposure^ in the 
Y, J, H and K bands. The resulting magnitude limits are typi- 
cally Fvega - 20.2, 7vega - 19.6, //vega - 18.8 and /Tvega - 18.2 

for point- sources detected wi th S/N^ 5, and the average see- 
mg IS 0'.'8 dWarren et al.ll2007l) . The area, depth and wavelength 
coverage of the LAS were chosen with the detection of z - 6 
quasars in mind, and the Y filter, positioned between the z and J 
bands in a region of minimal atmospheric absorptio n (0.97 //m % 
A % L07 jum), is particularly useful in this regard dHewett et al.l 
2006). Quasars in the redshift interval 6.4 ^ z ^ 7.2 would be 
identifiable as being unusually red in i-Y and z-Y, while be- 
ing significantly bluer in Y-J than the more numerous L and T 
dwarfs. 

All UKIDSS images and catalogues are made available to 
European Southern Observatory (ESO) countries and then, 18 
months later, the world, in a series of incremental data releases. 
The re sults presented here are based on the Third Data Release 
(DR3; IWarren et al.ll2009l) for which there is 870 deg^ of LAS 
coverage in the Y and J bands. 



3. Candidate selection 

High-redshift quasar candidates are selected in a two-step pro- 
cess, using survey data initially (Section 13.1b . followed by ad- 
ditional photometric observations of the most promising objects 
(Section[l2Ti. 

3. 1 . Initial shortlist 

Redshift ~ 6 quasars are expected to appear as stationary point- 
sources with extremely red optical-NIR colours. The first stage 
of the candidate selection process is to extract a fairly complete 
(if highly contaminated) sample of all such sources from the 
UKIDSS and SDSS databases. 

The initial selection from the UKIDSS LAS was of all point- 
sources with ivega ^ 19.88 that were also detected in the J 
band and have measured Fvega-./vega 0.88 (to exclude L and 
T dwarfs). The resulting sampl e was then cross-matched to the 
SDSS Fifth Data Release (DR5;'Adel man-McCarthv et al.l2007l) 
and all sources either undetected by SDSS or with iAB-^Vega ^ 
2.5 (but undetected in u, g and r) were selected. Any sources with 
inter-band (or inter-survey) positional mis-matches of greater 
than 0'.'7 were rejected to avoid nearby Galactic stars with ap- 
preciable proper motions, as well as most asteroid^. 

A number of heuristic algorithms were applied to reject 
sources for which the UKIDSS or SDSS database photometry is 
likely to be unreliable. Contaminating sources elimin ated at this 
stag e include those arising from WFCAM cross-talk dDve et al.l 
2006), objects in the haloes of bright stars, and close pairs sepa- 
rated by a few arcsec for which deblending is required in SDSS. 



^ The WS A is located at http : / /surveys . roe . ac . uk/wsa/. 

* 80 s exposures are sometimes taken in mediocre conditions, an ex- 
ample of which is the finding chart shown in Fig.[T] 

^ Some asteroids, observed at a turning point between prograde and 
retrograde motion, can appear stationary in UKIDSS survey images 
taken within an hour of each other 
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For the sources that were undetected in SDSS i and z, aper- 
ture fluxes, corrected for aperture loss, were measured from the 
SDSS images. Acquisition of the aperture fluxes represents the 
final step in a fully automatic procedure that, in the analysis of 
UKIDSS DR3, yielded -2x10"* pre-candidates with /-, z-, Y- 
and 7-band photometry. 

The next stage of the selection process is to use Bayesian 
model comparison to determine the probability, given the photo- 
metric data, that each source is a hig h-redshift qua s ar, P g . This 
calculation is described in detail in Mortl ock et a D (m)9). The 
key to this approach is having an accurate model of the stellar 
population, specifically the intrinsic magnitude and colour dis- 
tributions of the cool M dwarfs which scatter into the quasar 
selection box described above. Convolving this distribution in 
z, Y and / with the observational noise gives the likelihood 
that any candidate is an M dwarf. This can be compared directly 
with the likelihood that the candidate is a quasar, which is cal- 
culated using an extrapolation of the high-redshift quasar lumi- 
nosity function of Fan et al. (2004) and simulated colours from 
iHewett et alj (12006 ). The result is a fully self-consistent value 
of f q that, in principle, combines the available information on 
the candidate with the constraints on the quasar and star popula- 
tions in an optimal way. Of course there are practical limitations, 
the most important of which result from incomplete sampling of 
the tails of the noise distributions and the limited knowledge of 
the stellar population fainter than the UKIDSS and SDSS survey 
limits. The latter is important as the probability calculation nec- 
essarily includes the possibility of faint sources scattered into 
the sample, but is difficult to assess empirically using the sur- 
ve y data. Despite these ambiguities (which are explored further 
in lMortlock et al.|[2009!) , the probabilities which result from the 
adopted noise and population models clearly provide a good ob- 
jective ranking scheme for quasar candidates. 

Most of the candidates are clustered close to the stellar lo- 
cus and their low quasar probabilities (i.e., Pq $ 0.01) merely 
confirm the obvious. The motivation for the method is to assess 
the ambiguous sources lying between the quasar and star loci 
in colour space, or close to the F-band limit with quasar-like 
colours but large errors. The sometimes counter-intuitive results 
are explored in some detail in Mortlock et al. (2009), but the crit- 
ical point is that the candidates can be ranked by Pq and thus 
prioritised objectively for further investigation. 

It is only at this stage that any visual inspection was un- 
dertaken, with the UKIDSS and SDSS images of the few hun- 
dred best candidates checked for artefacts such as hot pix- 
els, bad columns, obviously wrong photometry, and undetected, 
blended or moving sources. (Although the UKIDSS and SDSS 
data reduction pipelines flag the vast majority of such in- 
stances, the selection of unusual sources on the basis of ob- 
served inter-survey colours inevitably produces samples with 
an over-representation of the very rare cases for which at least 
one of the surveys has produced anomalous measurements.) 
Following visual inspection, only a small sample of S 100 plau- 
sible quasar candidates with f q <: 0.01 remainecQ, one of which 
was ULAS J1319+0950, shown in Fig. [T] Follow-up observa- 
tions were needed to either reject these candidates as stars or, 
hopefully, confirm some as high-redshift quasars. 



* The choice of threshold value for Pq (which defines the list of candi- 
dates and determines the completeness and contamination of the quasar 
sample) depends on the resources available for follow-up photometry 
and spectroscopy. 




Fig.l. Finding chart for ULAS J131911.29-I-095051.4 showing 
the UKIDSS LAS 80s /-band image of a field 3 arcmin on a 
side centred on the source position at RA - 1 99.79706 deg and 
dec = 9.84762 deg. 



3.2. Follow-up photometry 

Rather than immedia tely taking spectra of all the candidates (cf 
iGlikman et al."2008'), it is more efficient first to refine the pho- 
tometry, obtaining short exposures in the /, z, Y and J bands 
on a variety of telescopes. The quasar probability defined in 
Section lTTl is recalculated whenever new data are obtained, and 
a source is dropped from the candidate list if Pq falls below the 
selection threshold at any point. As almost all the candidates ini- 
tially have minimal S /N in the / band SDSS images, follow-up 
observations in this band are the most efficient way to quickly 
reject such sources, revealing most to be cool stars scattered to 
have quasar-like colour by photometric noise. Observations in / 
are thus prioritised, although initial follow-up observations are 
sometimes made in z, Y or J, depending on weather, telescope 
scheduling, and other external factors. Candidates are not con- 
sidered for spectroscopy until follow-up photometry has been 
obtained in at least the /, Y and J bands. 

Follow-up observations of ULAS J13 19+0950 were ob- 
tained in /, z, Y and J, as summarised in Table [1] 

ULAS J1319-I-0950 was first re-observed in the i band at the 
Liverpool Telescope (LT) on the nights beginning 2008 January 
8 and 12 for a total of 3240 s. The candidate's quasar proba- 
bility remained high, and so 300 s exposures in each of the Y 
and J bands were obtained using the UKIRT Fast-Track Imager 
(UFTI) on the night beginning 2008 January 16. The candidate 
still appeared promising so the UFTI observations were repeated 
the following night. 

The improved /-, Y- and 7-band measurements were suf- 
ficiently precise that ULAS J1319+0950 had Pq ^ 1 (i.e., it 
was essentially implausible for a cool star to have scattered to 
the candidate's observed colours given the precision of the new 
photometric data), whereas all the other initial candidates were 
rejected on the basis of their follow-up photometry. A spec- 
ti-um of ULAS J1319+0950 was thus obtained on the night 
beginning 2008 January 22, with the result that it was imme- 
diately confirmed as a high-redshift quasar (Section |4]i. That 
the entire selection, follow-up and confirmation process for 
ULAS J1319+0950 was completed less than seven weeks after 
the UKIDSS DR3 release (on 2007 December 6) is a powerful 
illustration of the advantages of queue-scheduling for the three 
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telescopes involved, especially considering the RA of the target 
and the time of year. 

ULAS J1319+0950 had already been detected as a faint 
source in SDSS, although with zab = 20.13 ± 0.12 its S /N was 
too low for it to be selected in to the high-redshift quasar sam- 
ple defined bv lFan et al.l (|2003|) . Post-confirmation, improved z- 
band photometry of ULAS J1319+0950 was obtained using the 
ESO Multi-Mode Instrument (EMMI) on the New Technology 
Telescope (NTT) on the night of 2008 January 29. The observa- 
tions were made with the long -pass #611 filter, the bandpass of 
which, in combination with the red cut-of f of the CCD response , 
is quite similar to the SDSS z bandpass ( Vene mans et al.ll2007l) . 
To calibrate the NTT-image, SDSS /- and z-band photometry of 
bright, unsaturated stars in the frame was converted to the natu- 
ral system of the image using Zntt.ab = Zab - 0.05(/ab - Zab) 
dVenemans et al. 2 007) . and the result quoted in Table[T]is in this 
natural system. 

4. ULAS J1 31 9+0950 

After the series of follow-up photometric observations described 
in Section |T2l showed ULAS J1319-H0950 to be a promising 
high-redshift quasar candidate, an optical spectrum was obtained 
to confirm the identification (Section |4Tt . A more accurate red- 
shift was estimated from a NIR spectrum covering the Mgn 
emission line (Section 14.2b . after which both spectra were used 
to compare ULAS J1319+0950's emission and absorption prop- 
erties with those of similarly distant quasars discovered in SDSS 
(Sectionl431l. 

4.1. Spectroscopic observations 

An optical spectrum of ULAS J1319+0950 was obtained using 
the Gemini Multi-Object Spectrograph (GMOS) on the Gemini 
South Telescope on the night beginning 2008 January 22. Two 
spatially-offset spectra, each of duration 900 s, were obtained us- 
ing a r.'O slit and the R400 grating, covering the wavelength 
range of 0.5-1.0//m over the three CCDs. The standard bias sub- 
traction and flat-fielding steps were followed. Then, because of 
the strong sky lines in the red part of the spectrum, the 'double 
subtraction' method was used for the first-order sky subtraction 
(i.e., frame B was subtracted from frame A, and then the neg- 
ative spectrum subtracted from the positive spectrum, after the 
two spectra were aligned). This procedure removes most sys- 
tematic errors, but at the price of increasing the noise in the fi- 
nal frame by a factor of ~ 1 .4 compared to the theoretical limit. 
Second-order sky subtraction was achieved by fitting a smooth 
function to each column. At this point cosmic rays were removed 
using the Laplacian Cosmic Ray Removal Algorithm (LCRRA; 
Ivan Dokkumll200i1) . Wavelength-calibration was carried out us- 
ing observations of a Cu Ar lamp. Relative spectrophotometry, 
and correction for telluric absorption, was achieved using obser- 
vations of a standard star, and the spectrum was then scaled to 
match the NTT z-band photometry given in Table [T] The final 
GMOS spectrum is shown in Fig.|2] 

The spectrum is recognisable as that of a z - 6 quasar from 
the presence of a broad emission line, identified as Ly a, at the 
same wavelength as a strong continuum break, attributed to Ly a 
forest absorption. The Ly a emission peaks at a wavelength of 
A ^ 0.88 fim, which leads to a preliminary redshift estimate of 
z = 6.12. However the Ly a line is strongly absorbed to the blue 
and such high-ionization lines can ex hibit velocity shifts relative 
to the quasar systemic redshift (e.g., iTvtler & Fanlll992h . limit- 
ing the utility of the optical redshift measurement. 




0.6 0.7 0.8 0.9 

wovelength (/um) 

Fig. 2. The Gemini GMOS spectrum of ULAS Jl 3 19+0950 
(black curve) and the noise spectrum (red curve), both binned 
by a factor of four. The gaps at wavelengths of ~ 0.68 //m 
and ~ 0.82jum correspond to the breaks between the different 
CCDs. The wavelengths of common emission lines redshifted 
by z = 6.13 are indicated. 

The Mgn line, if observable in the K band, should give 
a more reliable estimate of the redshift, and so a NIR spec- 
trum of ULAS J1319+0950 was obtained. The source was ob- 
served using the Near-IR Instrument (NIRI) on the Gemini North 
Telescope on the two nights beginning 2008 February 26 and 
27. Observations were made with a 0'.'75 slit and the K grism 
G5204, covering the wavelength range 1 .9-2.5 fim with a resolv- 
ing power of R - 500. With NIRI it is standard procedure to dis- 
card the first exposure of a sequence, leaving a total of 1 1 usable 
300 s exposures over the two nights. The o bservation and data 
reduction methodologv of Weatherlev et al. (2005) was adopted, 
with six different slit positions used to reduce the noise from 
sky-subtraction relative to the more standard ABBA sequence. 
The data suffered from increasing slit losses as the observations 
proceeded, which is believed to be due to differential flexure be- 
tween the guide probe and the instrument. The additional slit 
losses, relative to the first frames on each night, reduced the final 
S /N of the detected emission line by a factor of ~ 1.5. Because 
of the varying throughput, each sky-subtracted frame was first 
scaled to a common count level, and then weighted by the in- 
verse variance in the sky (as evaluated in a region free of strong 
emission lines). Wavelength cal ibration was performed using the 
list of sky emission lines from iRousselot et al] (l2000h . A stan- 
dard star observed at similar airmass was used to correct for tel- 
luric absorption, and for relative flux calibration, after which the 
spectrum was scaled to match the UKIDSS A'-band photometry 
given in Table [1] The final NIR spectrum is plotted in Fig. [3] 

4.2. Redsliift estimation 

The NIR specti-um of ULAS J 1 3 1 9+0950 shown in Fig.|3]reveals 
the broad Mg ii emission line near /I ^ 2.0 fim, corroborating the 
initial redshift estimate of z - 6.1. To obtain a more accurate 
measurement (and to quantify the absorption in the Ly a forest 
in Section I4.3l l a power-law fit to the continuum was made to 
the combined optical and NIR spectrum by minimizing and 
iteratively clipping outliers to eliminate the emission lines. The 
best-fit power-law is /)(/l) oc /l"*^""*, where a - 0.32 + 0.01 ( de- 
fined so that fy oc v^"). This is plotted in Fig.|4] This continuum 
was subtracted from the data and then a Gaussian was fit to the 
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Table 1. Original survey and follow-up photometric observations of 
ULAS J1319+0950. 



filter 


original 


follow-up 


observing date 


telescope 


exposure time 


'AB 


22.83+0.32 


22.55+0.09 


2008 January 8 and 12 


LT 


2x 1620 s 


Zab 


20.13±0.12 


19.99+0.03 


2008 January 29 


NTT 


450 s 


IVega 


19.22±0.06 


19.10+0.03 


2008 January 16 and 17 


UKIRT 


2 X 300 s 




18.69±0.05 


18.76+0.03 


2008 January 16 and 17 


UKIRT 


2 X 300 s 



All magnitudes are quoted in the natural system of the initial survey: AB for the optical SDSS 
bands and Vega for the NIR UKIDSS bands. Note that the transmission profile of the zntt filter 
dilfers somewhat from that of the SDSS z filter (see Section [J!2l (. 



1.0 n 



0.8 - 




0.2 - 



0.0 

1.95 2.00 2.05 2.10 2.15 2.20 
wavelength (/xm) 

Fig. 3. The Gemini NIRI spectrum of ULAS J1319H-0950 (black 
curve) and the error spectrum (red curve), both binned by a 
factor of two. The continuum -t-line fit to the Mgii emission 
line is also shown (blue). The fit implies a source redshift of 
z - 6.127 ± 0.004. The effect of telluric absorption was removed 
using a standard star observed at similar airmass, resulting in the 
increased errors at ~ 2.005 yum. 

residual Mg ii emission line; the resultant continuum+line fit is 
plotted in Fig. [3] The central wavelength of the Mg ii emission is 
A = (1.994 + 0.001)yum, implying that ULAS J1319-I-0950 has 
a redshift of z = 6.127 + 0.004. The Mgii emission line has a 
rest-frame equivalent width of ~ 14 A, substantially less than the 
typica l value of ~ 28 A for z - 6 quasars found by iKurk et al] 
(l2007h . 

4.3. Comparison with SDSS quasars 

ULAS J1319+ 0950 is, after ULAS J0203 -1-0012 
jVenemans et al.l 12007). only the second redshift ~ 6 quasar 
discovered in the NIR. Given that surveys at these wavelengths 
are the first capable of probing beyond the current limit of 
z - 6.4, it is important to determine whether quasars selected 
in the NIR exhibit significant differences from those selected 
in the optical. This preliminary exploration focuses on the 
emission line strengths, the neutral hydrogen optical depth, and 
the ionization region around the quasar 

4.3.1. Emission line properties 

Redshift ~ 6 quasars are identifiable in photometric surveys 
primarily due to their extreme colours which result from the 
strong neutral hydrogen absorption blueward of the quasars' 
Ly a emission lines. However Fig. |5] which compares the rest- 
frame spectrum of ULAS J1319+0950(and ULAS J0203+0012; 




0.5 1.0 1.5 2.0 2.5 

wavelength (/tim) 



Fig. 4. The optical spectrum, binned by a factor of four (left 
black curve), and the NIR spectrum, binned by a factor of five 
(right black curve), of ULAS J1319-H0950. The power-law con- 
tinuum fit is also shown (blue). 

see Section|5]) to a composite SDSS quasar spectrum ("Fa n et all 
[2006b), shows that ULAS J1319+0950 has a noticeably weaker 
Ly a emission line than is typical for optically-identified sources. 

Measuring the rest-frame equivalent width of 
ULAS J1319-h0950's Lya line is somewhat problematic 
as the continuum level is ambiguous and the power-law fit 
shown in Fig. |4] is probably too steep to give a good local 
continuum estimate. For the purposes of comparison the contin- 
uum levels shown in Fig. |5] were adopted, yielding rest-frame 
equival ent widths of 19 3 A for ULAS J1319+0950 and 60.2A 
for the iFan et al.1 (|2004|) composite. This quantifies the visual 
impression that ULAS J1319-I-0950 has a weaker Ly a emission 
line than an average z - 6 quasar, although there are several 
examples of h igh-redshift quasars with even weaker lines (e.g., 
iFan et al.|[T999; Dia mond -Stani c et al . 2009b . 

If ULAS J1319-H0950 differed only in that its Lya emis- 
sion was marginally stronger (or, equivalently, less absorbed to 
the blue), it would almost certainly have been discovered in the 
SDSS high-redshift quasar sample. With zab = 20.13 ± 0.12 
in the SDSS database, it was brighter than the iFan et al] (|2003|) 
limit of Zab - 20.2, but it failed the S/N cut of cr^^g < 0.1. 
Given that the NTT re-observation described in Section [J!2] gave 
Zab = 19.99 + 0.03, the SDSS team were unlucky that the mea- 
sured z-band flux was less than the true value and that the source 
was observed in slightly worse than average conditions, both of 
which combined to leave the source outside their selection cuts. 
This illustrates the more general limitation inherent in any selec- 
tion method which applies a hard data cut in a region of parame- 
ter space where there are appreciable numbers of objects. A way 
to avoid the problem, at least in principle, is to apply the proba- 
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0.15 



Fig. 5. Rest-frame spectra of ULAS J0203+0012 (blue) and 
UL AS J 1 3 1 9+0950 (green) compare d to the composit e spectrum 
of z ^ 6 SDSS quasars presented by lFan et al ] (12004 . The con- 
tinuum level adopted for the calculation of the equivalent with 
of the Ly a line for each spectrum is also shown (dashed lines). 

bilistic approac h adopted in Section [TTI and described in full by 
iMortlock et al.l (12009.) . although it remains difficult to deal with 
the large numbers of low S /N candidates near the survey Umit. 

4.3.2. Neutral hydrogen optical depth 

Observations of z ^ 6 quasars have revealed a marked increase 
in th e density of ne utral hydrogen above a redshift of ~ 5.7 
(e.g.. lFan et al.ll2002l), possibly indicating the end of cosmologi- 
cal reionization (e.g.. lFan et al.ll2006a) . As hydrogen absorption 
also has an effect on the detectabiUty of z - 6 quasars, it is im- 
portant to assess whether the measured optical depths of optical- 
and NIR-selected quasars differ significantly. 

Here, the spectru m of ULAS J13 19+0950 is analysed fol- 
lowing the method of iFan et alj (l2006at) . with the effective opti- 
cal depth, TefF(zabs), St redshlft Zabs estimated to be 



Teff (Zabs) = - In 



</<,obs[(l +z)^Lya]);^;,i,, 
</),mt[(l +z)^Lya]>z-z,b., 



(1) 



where Xi,obs(^) is the observed flux density, /),int('t) is the con- 
tinuum flux density (assumed to be given by the power-law fit 
described in Section |4!2] ). and Ai^a = 0.1216/zm is the rest- 
frame wavelength of the Ly a transition. The angle brackets de- 
note an averaging over a finite redshift range around Zabsi bins 
of width Azabs = 0.15 were used. Given that the Lyy6 and Lya 
emission lines appear at A - 0.74jum and A = 0.87 yum, respec- 
tively, Teff(Zabs) Can be estimated in the range 5.1 ^ Zabs ~ 6.0. 
The errors in these optical depth estimates are dominated by the 
noise in the spectrum of ULAS J1319+0950, with the uncer- 
tainties in the continuum fit and sample variance (across differ- 
ent lines-of-sight) being secondary. The resultant error bars are 
also significantly asymmetric, especially for the Zabs = 5.93 bin, 
in which flux is detected with only S/N^ 2. These Ly a optical 
depth estimates, along with an analogous estimate of the Ly/3 
op tical depth (with the Ly a emission corrected for as described 
bv lFan et al.ir2006al) . are given in Table |2] 

These estimated optical depths towards ULAS Jl 3 19+0950 
are c onsistent with measurements of SDSS quasars by Fan et alj 
2006a), which gives some confidence that Gunn & Peterso^ 
1965h effect measurements of optical- and NIR-selected z - 6 



Table 2. Effective optical 
depths of the absorption sys- 
tems seen in the spectrum of 
ULAS J1319+0950. 



transition Zabs 



Teff 



Lya 
Lya 
Lyor 
Lya 
Lyor 
Lya 

LyyS 



5.93 
5.78 
5.63 
5.48 
5.33 
5.18 



T, «:q+0.7l 
A 97+0.45 

T TC+0.13 



2.58+ 
3.21! 



^■^^-0.03 



7 4^+0.06 



The Ly /? optical depth has been 
corrected for Ly a absorption 
usi ng the rela tion found by 
iFan et aLll [2006al) . 



quasars may be combined. 



4.4. Quasar ionization region 

The ultraviolet (UV) photons from quasars are sufficiently ener- 
getic to i onize the neutral h ydrogen in a volume around them, 
creating a'Stromgren] (Il939h sphere. The extent of the ionization 
region depends in part on the UV luminosity of the quasar and 
the neutral fraction of the surrounding medium, potentially pro- 
viding a complementary probe of reionization. Unfortunately, 
the measurement is difficult in practice, due to uncertainties in 
the quasars' UV spectral energy distributions (SEDs), unknown 
ionizing flux from nearby galaxies, and difficulty in establishing 
the extent of the H ii region from hea vily absorbed and otherwise 
noisy spectra (e.g.. lFan et al.ll2006al) . 

It has proved more useful to adopt a relative approach in 
which the evolution of the size of quasars' Hii regions with 
redshift is estimated without reference to an absolute model of 
the neutral hydrogen fraction. The relative measurement is at- 
tempted here fo r ULAS J 1319+0950, once again following the 
methodology of lFan et akl (j2006a), by defining the quasar prox- 
imity region as the volume within which the transmitted flux 
fraction is $ 0.1. After dividing the Gemini spectrum by the 
continuum+line fit (described in Section |4~2l ) to convert to trans- 
mitted flux, the spectrum was smoothed to a resolution of 20 A. 
This smoothed transmission spectrum remains below 0. 1 down 
to a wavelength of A ^ 0.856/im; converting to redshift and then 
to a co-moving physical length implies that ULAS J1319+0950 
has a proximity zone of radius ^ 5.1 Mpc. 

The inferred radius is broadly consiste nt with m easurements 
from optically-detected SDSS quasars (.Fan et al.l [2006a). al- 
though to make a more quantitative comparison it is neces- 
sary to standardise the above measurement of Rj, to a fidu- 
cial absolute magnitude Mi45()_ab = -27 according to Rp^27 = 

100-4(Ml450,AB+27)^p^ With Mi45().AB = -27.12 this gives Rp^27 ^ 

4.9 Mpc for ULAS J1319 +0950. This is low er than the fiducial 
value of ~ 7 Mpc given bv lFanetalJ(l2006al) forz - 6.1 quasars, 
although this difference is not significant given the range of Rp 2i 
values seen in the optical sample. 



5. ULAS J0203+0012 

ULAS J0203+0012wa s die first z ^ 6 quasar found in UKIDSS, 
IVenemans et al.l ( l2007h describing the discovery and reporting a 
re dshift of z = 5.86 . The source was subsequently recovered 
by iJiang et all (l2008l) in their survey for high-redshift quasars 
using deep co-added data from multiple scans of SDSS Stripe 
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Fig. 6. Re-calibrated VLT spectrum of ULAS J0203+0012 
(black curve) and the noise spectrum (red curve), binned by a 
factor of two. The wavelengths of common emission lines red- 
shifted by z = 5.72 are indicated. 



82; they quote a redshift of z = 5.85. An improved redshift 
for ULAS J0203H-0012 was sought by acquiring a NIR spec- 
trum (Section imi. These data revealed ULAS J0203-H0012 to 
be a broad-absorption line (BAL) quasar, resulting in a modi- 
fied redshift estimate (Section |S!2l ), and also showing some dis- 
crepancies compared to the existing optical spectrum in the re- 
gion 0.9-1.0;um. The origin of the discrepancy was traced to the 
relatively poor quality of the published standard star data used 
to flux-calibrate the optical spectrum, which has hence been re- 
reduced (Section lSTt . 



5.1. Spectroscopic observations 
5.1.1. Optical spectrum 

The original optical spectroscopic observations of 
ULAS J0203-I-0012 were made in 2006 September using 
the FOcal Reducer and low dispersion Spectrograph (FORS2) 
on the Very Large Tele s cope (VLT), and are described in 
detail in IVenemans et al.l (l2007l) . The re-reduction followed a 
similar sequence to that described there, up to the point of flux 
calibration. Instead of using the published spectropho tometry of 
the white dwarf standard star GD50 from a model 

spectrum with effective temperature ^ eff - 41 150 K and su rface 
gravity log[g/(l cms"^)] = 9.15 from lDobbie et alj (l2005l) was 
used to establish the flux calibration curve. Finally, the data 
were scaled to match the measured z-band photometry. The 
revised spectrum, plotted in Fig. |6] shows a sharp downturn 
longward of 0.95 fim compared to the original spectrum. 



5.1.2. NIR spectrum 

A NIR spectrum of ULAS J0203H-0012 was obtained using the 
Gemini Near Infra-Red Spectrograph (GNIRS) on the Gemini 
South Telescope on the night beginning 2007 March 22. The 
observations were taken in cross-dispersed mode with the short 
camera, the 32 lines mm"' grism, and a 0'.'75 slit, providing 
coverage from ~ 0.8yum to ~ 2.5 yum at a resolving power of 
R - 500. Due to the limited slit length it was impractical to 
adopt the observation strategy used for the NIRI observations of 
ULAS J1319H-0950 (SectionO and so a standard ABBA off^set 
pattern was used. 




1.5 

wavelength (Aim) 
Fig. 7. The combined optical and NIR spectrum, spliced at /I = 
0.96 Aim, of ULAS J0203 -1-0012 (black curve) and the noise 
spectrum (red curve), both binned by a factor of eight. The 
power-law continuum fit is also shown (blue). The wavelengths 
of common emission lines redshifted by z = 5.72 are indicated. 
The Mgii line would be expected to appear at 1.88A'm, in the 
middle of the strong telluric absorption band at ~ l.9iJ.m for 
which the data have been omitted. 



The observations comprise eight frames, with a total expo- 
sure time of 2400 s. The data were reduced mainly using the 
Gemini-GNIRS package in the Image Reduction and Analysis 
Facility (IRAF; lTodvl[l993b . After correction for pattern noise 
and flat-fielding, each of the diffraction orders was separated, 
and the double-subtraction method (described in Section 14.11) 
was applied to each of the four pairs of images. The resulting 
four frames for each order were then averaged. S-distortion cor- 
rection and wavelength calibration were then applied using ob- 
servations of an Ar lamp, and one-dimensional spectra were ex- 
tracted. The data were corrected for telluric absorption and flux- 
calibrated using the spectrum of a spectroscopic standard star 
Finally, the different diffraction orders were spliced together to 
produce a single continuous spectrum. Strong telluric absorption 
bands occur at 1.35-L43Aim and 1.80-1. 95 A^m, and so the data 
in these wavelength ranges were discarded. The spectrum was 
placed on an absolute flux scale by matching to the /T-band pho- 
tometry. 

Comparison of the optical and NIR spectra in the region of 
overlap (Q.^-l.Qfim) showed good detailed agreement, except 
for a 1 per cent normalisation offset, and so the optical spectrum 
was scaled to the NIR spectrum. The declining S /N of the NIR 
spectrum towards the blue, and of the optical spectrum towards 
the red, are approximately equal at 0.96 A^m, so the two spectra 
were spliced at this point. The final combined spectrum is plotted 
in Fig.|7] 

5.2. Redshift estimation 

The most striking features of the combined ULAS J0203 -1-0012 
spectrum shown in Fig. [T] are the broad absorption lines near 
0.9 A^m and 1.0/vm. These are attributed to Siiv and Civ, re- 
spectively. Despite the increased wavelength coverage, the red- 
shift is difficul t to determine accu ra tely, as is common with BAL 
quasars (e.g. Trump et al.l |2006|) . IVenemans et al.l (12007 ) and 
Jiang et al. (2008) identified the peak of the line near 0.835 A^m 
with Ly a, and the shoulder at 0.855 A^m with N v, but matching 
Si IV, C IV and C iii] lines for this redshift are not evident. Instead 
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it seems likely that the peak at 0.835 //m is the N v emission line, 
and that the Ly a emission line has been absorbed by a N v BAL. 
As shown in Fig. |7] a redshift of z = 5.72 provides a reasonable 
match to the N v peak, as well as weak features that match Si iv. 
Civ, and Cm]; but any redshift in the range 5.70 < z < 5.74 
is consistent with the data. Unfortunately the Mg ii line, which 
ought to give a definitive redshift, lies in the region of atmo- 
spheric absorption near 1 .9 yum. 

In the wavelength region 0.85-1. 08 /^m it is difficult to de- 
termine the extent of the BAL troughs, and the regions of un- 
absorbed continuum emission. The best fit power-law, plotted 
in Fig. |7] is a poor match. It seems likely that the sharp step 
near 0.865 yum marks the blue edge of the Siiv BAL, and that 
the shoulder at 0.855yum, previously identified with Nv, is in 
fact continuum emission. The resulting estimate of the rest frame 
equivalent width of the N v line is ~4.2 A. 

Given the presence of BAL troughs, it is difficult to infer 
very much about the neutral hydrogen along the line-of-sight to 
ULAS J0203H-0012, as the Ly a emission line is all but com- 
pletely absorbed. Similarly, it is not possible to measure the 
iGunn & PetersonI ^1965) optical depth to ULAS J0203-H0012: 
assuming the N v BAL has a similar velocity width to the C iv 
and Si rv BALs, the absorption blueward of Ly a is most likely 
dominated by the wing of the N v BAL, and not by intervening 
neutral hydrogen along the line-of-sight. 

The case of ULAS J0203H-0012 is rather similar to that 
of the BAL quasar SDSS J1044-0125 jFane t al. 200C |. The 
initial redshift measurement of z = 5.80 by Fan et al.l (l2000t) 
was revised to z = 5.74 bv [Goodrich et al. ( 200 lb o n the ba- 
sis of NIR spectroscopy, after which lJiang etal. ' ('2007') used the 
Cm] line to estimate z = 5.78. Both ULAS J0203+0012 and 
SDSS J1044-0125 demonstrate the importance of NIR spectra 
to obtaining reliable redshifts of z ^ 6 quasars (and especially 
BALs). 



6. Conclusions and future prospects 

Two new high-redshift quasars have now been discovered in 
the UKIDSS LAS: ULAS J020 3-H0012 at z ^ 5.72 (pre- 
viously given as z = 5.86; IVenemans et al.1 l2007h ; and 
ULAS J1319H-0950 at z = 6.13. Two previou sly known high- 
redshift quasars, SDSS 0836H-0054 at z = 5.82 ("Fan et al.ll200l 
and SDSS 1411 + 1217 at z = 5.93 (Fan et al. 2004), have also 
been recovered in the UKIDSS DR3 LAS area. The two new 
quasars are well within the redshift range covered by optical sur- 
veys, but they were both too faint in SDS S to be include d in the 
high-redshift quasar sample defined by F an et al.l (l2003l) . While 
the revised lower redshift of ULAS J0203+0012 puts its Lya 
emission near the peak of the SDSS z band response, the line is 
almost completely absorbed by a N v BAL, which results in a re- 
duction of the z-band flux. ULAS J1319-I-0950 has a very weak 
Ly a line, with a rest-frame equivalent width of just ~20 A. The 
small Lya equivalent width results in the z-band S /N falling 
ma rginally b e low th e limit of the SDSS quasar sample defined 
by iFan et al.] (l2003l) . Nonetheless, the detection of four z - 6 
quasars brighter flian Fvega = 19.88 in the 870 deg^ of UKIDSS 
DR3 LAS is consisten t with expectations from the luminosity 
function of iJiang et al.l (12008) . and so these new detections do 
not imply the existence of a previously unidentified quasar pop- 
ulation with weak emission features. 

These results demonstrate that the UKIDSS data are of suffi- 
cient quality to recover z ^ 6 quasars down to F ^ 19.5 with rea- 
sonable completeness. Quantification of the completeness will 



be explored in more detail when the sample is larg er. At the time 
of writing the UKIDSS Fifth Data Release (DR5; 'Warr en etal.l 
2009) has already taken place, increasing the LAS area covered 
in the Y and / bands to ~ 1360deg~. DR5 is of sufficient size 
that finding no quasars of redshift z > 6.4 would be somewhat 
inconsistent with expectations. 
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